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INTRODUCTION  WITH  GOALS 


The  goal  of  our  segment  of  the  Unified  Physical  Parameterization  for  Extended  Forecasts  DRI 
continues  from  prior  reporting  periods  to  provide  observations  related  to  the  partitioning  of 
atmospheric  boundary  layer  (ABL)  fluxes  into  two  contributions:  eddy  diffiisivity  and  vertical  mass 
fluxes  (EDMF).  In  addition  to  the  field  campaign,  the  goal  is  to  validate  the  existing  EDMF 
expressions  and  to  propose  alternatives  or  modifications  to  those  formalisms. 

These  goals  are  being  pursued  using  data  sets  obtained  from  the  CIRPAS  Twin  Otter  flown  in  two 
differing  regions:  off  the  coast  of  California  near  Monterey  and  over  portions  of  the  Dugway  Proving 
grounds  in  Utah.  Data  provided  by  UCI  (Khelif)  from  the  Controlled  Towed  Platform  (CTV)  has  also 
been  incorporated  into  our  investigations. 

OBJECTIVES 

The  general  objective  is  to  collaborate  with  NRL  (Wang),  UCI  (Khelif)  and  the  general  numerical 
modeling  community  working  with  the  WRF  and  CO  AMPS  to  use  airborne  sensors  to  probe  the  ABL 
for  organizing  structures  such  as  OLEs  and  LLJs  and,  using  models  when  appropriate,  construct 
algorithms  that  would  relate  space-based  sensor  observables  to  the  varying  structure  dependent  ABL 
fluxes.  Two  regions  have  been  used  to  collect  data  for  the  evaluation  and  modification  of  the  EDMF 
formalisms: 

1.  In  areas  of  neutral  to  stable  ABLs  capped  by  cloud  cover  ranging  from  zero  to  100%,  measure 
the  fluxes  at  the  air-sea  interface  using  the  TODWL  and  CTV  instruments  flown  together  near 
the  California  coastline  in  September  2012. 

2.  Using  the  MATERHORN  experiment  at  Dugway,  UT  in  October  2012  as  an  opportunity  to 
expand  the  DWL  data  base  to  include  ABL  surveys  in  complex  terrain  as  a  comparison  to  the 
ABL  over  cold  water.  Note  that  the  MATERHO^  was  co-funded  by  ONR  (Ferek). 

As  this  research  has  progressed,  we  have  added  a  modeling  component  (See  attachment  4)  to  better 
understand  the  implications  of  ABL  structures  being  revealed  by  the  TODWL  during  the  September 
and  October  2012  field  campaigns. 

GENERAL  APPROACH 

The  work  on  objective  one  above  continues  as  we  examine  several  cases  in  2012  where  the  CTV  was 
operating  adequately  and  the  ABL  contained  OLEs  and  LLJs  needed  for  our  investigation.  The  first 
important  step  was  to  navigate  the  data  from  the  two  systems  so  that  their  data  sets  could  be  joined. 
This  was  not  as  easy  as  it  seemed  since  the  CTV  did  not  fly  faithfully  directly  below  the  Twin  Otter. 
This  platform  relative  navigation  has  been  done  and  is  reported  in  attachment  3. 

The  second  step  is  to  compute  heat,  moisture  and  momentum  fluxes  for  both  the  CTV  and  Twin  Otter 
flux  sensors.  The  computation  of  fluxes  from  the  CTV  will  be  done  by  both  UCI  (Khelif)  and  SWA 
(Emmitt). 


The  third  step  is  to  determine  the  correlations  between  ABL  structures  and  the  fluxes.  The  end  product 
is  envisioned  as  being  a  set  of  PDFs  for  mass  flux  contributions  as  a  function  of  regional  model 
representations  of  the  BL  including  space-based  observations  that  can  be  used  to  infer  organized 
structures  such  as  LLJs  (ASCAT  and  CALIPSO)  and  OLEs  (MODIS,  VIIRS  and  GEOS). 

The  fourth  step  is  to  evaluate  several  modifications  to  the  EDMF  as  currently  being  implemented  in 
models.  These  modifications  would  attempt  to  improve  the  representation  of  boundary  layer  roll 
(convective  and  dynamically  driven)  contribution  to  the  total  energy  and  mass  fluxes. 

A  final  effort  will  be  made  to  relate  the  fluxes  measured  at  flux  tower  sites  with  those  measured 
onboard  the  Twin  Otter  during  flybys.  During  those  passes  the  winds  and  aerosols  were  mapped  by  the 
TODWL.  These  data  sets  will  be  explored  for  relationships  between  the  BL  organization  and  flux 
intensities. 

WORK  COMPLETED 

In  the  final  half  of  year  3/4  on  this  project,  we  have  processed  more  than  50  hours  of  TODWL  data 
from  both  the  Monterey  and  DPG  areas.  We  have  combined  the  information  from  the  TODWL,  CTV 
and  Twin  Otter  sensors  to  establish  the  relationship  between  the  local  fluxes  and  the  energetic  of  LLJs 
and  OLEs.  The  WRF  model  has  been  setup  and  run  for  the  7  MATERHORN  missions  and 
comparisons  between  the  TODWL  wind  profiles  and  the  model  profiles  have  been  completed  for  one 
day. 

The  Triad  Interaction  Model  for  OLEs  has  been  modified  by  Foster  to  investigate  “stacked  OLEs”  seen 
in  the  April  2007  datasets  (see  reported  results  in  attachment  6) 

We  continue  to  modify  the  EDMF  (Eddy  Difflisivity  and  Mass  Flux)  parameterization  to  account  for 
the  differences  between  thermally  driven  convection  and  dynamically  driven  vertical  transports. 

Emmitt  and  de  Wekker  are  using  the  WRF  and  COAMPs  models  to  test  out  sensitivities  to  changes  in 
the  EDMF  related  to  our  field  data.  We  have  received  several  versions  of  the  EDMF  from  Joao 
Teixeira  for  testing. 


RESULTS 

Most  of  the  results  of  our  last  year’s  research  effort  were  reported  at  the  DRI  Workshop  in  Monterey  in 
August  2014,  the  annual  AGU  meeting  in  December  2013  and  at  a  SPIE  Conference  in  September 
2014  in  Amsterdam.  Five  of  these  presentations  are  appended  to  this  report: 


Boundary  Layer  Rolls  Observed  Above  and  Below  a  Jet  in  a  Marine  Boundary  Layer;  presented  at 
AGU  2013  by  Ralph  C.  Foster;  George  D.  Emmitt;  Kevin  Godwin;  and  Steven  Greco 

Investigating  the  impacts  of  LLJs  and  OLEs  on  ABL  exchanges  and  transports  using  an  airborne 
Doppler  wind  lidar;  presented  at  AGU  2013  by  G.D.  Emmitt ,  R.C.  Foster^,  K.  Godwin^  and  S.  Greco* 


Vertical  transport  of  TKE  within  Organized  Large  Eddies;  presented  at  the  DRI  workshop  2014  by  G. 
D.  Emmitt  (SWA),  S.  de  Wekker  (UVA)  and  R.  Foster  (Uwash) 


EDMF  in  the  Near-Neutrally  Stratified  Boundary  Layer;  presented  at  the  DRI  workshop  2014  by  R. 
Foster  and  G.  D.  Emmitt 

Investigation  of  coastal  zone  complex  boundary  circulations  and  their  potential  impact  on  energy 
transport;  presented  at  SPIE  Defense  and  Security  Conference  in  Amsterdam  2014  by  G.  D.  Emmitt, 
K.  Godwin  and  R.  Foster. 


In  addition  to  the  presentations  provided  as  appendices,  individual  reports  by  Ralph  Foster  and  Stephan 
de  Wekker  are  also  attached.  The  material  presented  by  de  Wekker  combines  MATERHORN  research 
results  that  were  funded  by  ONR,  ARO  and  NSF.  This  work  lays  the  foundation  for  UPP  related 
investigations  that  will  be  carried  out  in  Year  4. 

RELATED  PROJECTS 

ONR  contract  to  study  the  utilization  of  Doppler  wind  lidar  (DWL)  data  to  quantify  the  contribution  of 
organized  large  eddies  (OLEs)  to  fluxes  in  the  marine  boundary  layer  (DYNAMO  DRI  (S.  Harper, 
TPOC).  This  funded  effort  is  now  complete. 

PLANS  FOR  YEAR  5 

We  have  three  primary  goals  for  year  5  of  this  effort. 

1 .  Finish  the  application  of  a  methodology  (developed  last  year)  for  joining  airborne  DWL  data 
and  in  situ  flux  data  to  achieve  an  understanding  of  how  the  organized  ABL  structures 
modulate  ABL  fluxes. 

2.  Develop  an  expression  for  mass  fluxes  by  OLEs  in  the  MBL  that  is  suitable  for  use  in 
numerical  model  parameterizations.  This  will  be  derived  from  flux  data  measured  directly  by 
the  CTV  (Monterey)  and  towers  (Dugway)  and  associated  with  organized  ABL  structures 
mapped  from  the  DWL. 

3.  Use  the  COAMPS,  WRF  and  SCM  to  investigate  the  sensitivity  of  these  models  to  the  ranges 
of  values  for  the  “MF”  term  as  inferred  from  the  Monterey  and  Dugway  data  sets. 

We  will  pursue  publishing  our  results  should  the  findings  of  this  research  merit  this  effort. 


Attachment  1 


Boundary  Layer  Rolls  Observed  Above  and  Below  a  Jet  in  a  Marine  Boundary  Layer 

presented  at  AGU  2013 

by  Ralph  C.  Foster;  George  D.  Emmitt;  Kevin  Godwin;  and  Steven  Greco 
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I  Ralph  C  Foster^  G.  David  Emmitt^,  Kevin  Godwin^,  Steven  Greco^ 
J^University  of  Washington,  Seattle,  WA,  United  States. 
pSimpson  Weather  Associates,  Charlottesville,  VA,  United  States. 


•  Improve  physical  parameterization  of  Atmospheric  PBL  processes 

•  Local  vs.  non-local  fluxes  of  momentum,  heat  and  water  vapor 

•  Convection 

•  highly-skewed,  very  unstable  stratification 

•  Organized  Large  Eddies  (OLE  or  rolls) 

•  Less  Skewed,  near-neutral  to  moderately  unstable  stratification 

•  Aircraft  and  towed  platform  observations 

•  Doppler  Wind  Lidar  to  directly  measure  PBL  structure 

•  Flight-level  and  near-sutface  fluxes 

•  See  also  Poster 

•  Investigating  the  impacts  of  LLJs  and  OLEs  on  ABL  exchanges  and  transports 
using  an  airborne  Doppler  wind  lidar  (A43A-0227) 


I  Research  supported  by  ONR  Marine  Meteorology 
(unified  Parameterization  Project  DRI  (Dr.  Ron  Ferek) 


Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 
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Mean  wind  at  flight  level  is  genera^y  from  the  north 


Case  Stutly ' 

•  Off  Monterey  Bay,  30  Sep,  2012 

•  Prospecting  DWL  mode 

■  -2°  forward  stare 

-  ~6km  surface  intercept  from  300  m  flight  level 
'  Line-of-Sight  wind  speed  (VLOS) 

•  Stacked  VLOS  ~2m  vertical,  50  m  horizontal  resolution 

•  ~0.1  m/s  VLOS  precision 

■  ~Cross-mean-wind  flight  legs 

-  Mean  state  vertical  profiles  from  Twin  Otter  ascent/descents 

•  Fluxes:  see  poster  A43A-0227 

Oh* * «.  J  vat  fufi.s  .Shown  S.tackod  OI  L  ": 

•  Multiple  cross-wind  legs  in  same  region 

•  Organized  Large  Eddies  found  throughout 

•  Mid-PBL  wind  maximum  (low-level  jet) 

•  Baroclinic  mean  profile  (land  PBL  advected  over  the  bay?) 

•  OLE  above  and  below  LLJ 

•  Sometimes  connected,  sometimes  not 

•  What  dynamics  causes  this? 
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•  Simplify  Problem: 

•  Neutral  stratification 

•  Shear  effects  only 

•  Omit  effect  of  stratification  above  jet 

•  No  surface  buoyancy  flux 

•  Elevated,  thin  baroclinic  layer 

•  Variable  K(z)  (akin  to  MRF/YSU) 

•  Nonlinear  stability  analysis 

•  Floquet  triad  interaction 
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2nd-ary  mode  in  above-jet  layer 
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•  Co-linear 

•  Based  on  dominant  mode 

•  Time-dependent  solution 
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•  A  PBL  jet  can  induce  "stacked"  OLE 

•  Co-linear  OLE  triads  generate  complex  PBL  structures 

•  There  is  a  weaker,  larger-scale,  oblique  above-jet  mode 

•  OLE  rolls  are  very  common 

•  Near-neutra! 

•  Low  Skewness 

•  Likely  significant  flux  contribution 


Attachment  2 

Investigating  the  impacts  of  LLJs  and  OLEs  on  ABL  exchanges  and  transports  using  an  airborne 

Doppler  wind  lidar 

presented  at  AGU  2013 

by  G.D.  Emmitt\  R.C.  Foster^,  K.  Godwin^  and  S.  Greco^ 


Investigating  the  impacts  of  LLJs  and  OLEs  on  ABL  exchanges 
and  transports  using  an  airborne  Doppler  wind  lidar 

G.D.  Emmitf,  R.C.  Foster^,  K.  Godwin^  and  S.  Greco"' 
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^-Simpson  Weather  Associates,  Charlottesville,  VA,  USA;  ^-APL,  Univ.  of  Washington,  Seattle,  WA,  USA;  ^-KSG  Science,  Starkville,  MS,  USA 


Objectives 

•  utilize  the  Twin  Otter  Doppler  Wind  Lidar  (TODWL) 
and  Controlled  Towed  Vehicle  (CTV)  to  take  ABL  wind 
and  flux  measurements  over  marine  and  coastal 
environments 

•  Extend  prior  investigations  {2001-2008)of  LUs  and 
OLEs  in  the  MBL  and  PBLs. 

•  Investigate  and  characterize  the  presence  of  rolls 
(OLEs)  at  the  boundaries  of  stratocum ulus  topped 
MBLs. 


•  Study  the  potential  impact  on  the  development  and 
implementation  of  the  EDMF  into  forecast  models. 


TODWL  Data  Description 


Attribute 

Prrfomranee 

Mrtrir 

Comm^fita 

LOS  rotolution  (ipplics  to 
vertioal  ptofilcs  of  3D  winds 
n  well) 

SOm 

Rnnge  resolution  to  Iwrd  targets  (ground  or 
dense  cloud)  can  be  belter  than  10  meters. 

U.V.W  resolution 

<  l0cin/« 

<  3ein/t  for  ttntionary  ground  bMcd 

Mwimuai  range 

6 -30  km 

Very  ^qscndenl  upon  lerosob 

Time  to  complete  full  step 
lUne  oonicel  seen  for  wind 

-20  6ec 

12  point  step  sUre  with  .5  -2  second  dwells 

Sampling  frequency 

100  Hz 

Integration  of  several  ihols  is  lypioal  to 
improve  range  performanoe 

Structure  Prospecting  with  TODWL 


•  Feature  prospecting  uses  a  very  shallow  angle  below 
the  horizon  (~  -1  -3  degrees  for  a  300m  flight 
altitude). 

•  Results  in  ~  2  m  vertical  resolution  and  50  m 
horizontal  resolution  with  ~10  meter  sliding  sample. 

•  It  takes  ~  40  seconds  to  profile  100  meters  below  the 
aircraft. 


aerosol/wind  structures  below  the  Twin  Otter 

•  Processed  Twin  Otter  instrument  and  CTV  data  for 
time  series  of  u,  v,  w,  q,  and  0. 

•  Match  up  times  and  their  features  from  the  TODWL 
and  CABIN  data  sets  near  flight  level  and  CTV  at  CTV 
cruise  levels. 


ClRPASTwin  Otter  with  CTV  below 
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Twin  Otter 
Probe  Data 


TODWL  LOS,  SNR  and  Spectral  Width 
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Towed  Vehicle 
Flux  Data 


Twin  Otter:  30-Sep-2012  tfi20:39  to  10:1S.69  UTC 
360| - ■ - ' - 1 - 


4  6  8  10  12  14 
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Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Summary  of  Segment  Statistics 


Flight 

Segment 

TODWL 

Ahitode 

CTV 

Altitude 

Heading 

TKE 

Sensible 

He8t(W) 

latent 

Heat(W) 

Skewness 

1007 

284 

60 

94 

.19 

8.S5 

3.43 

-.40 

1.92 

•1.28 

15.15 

.90 

1028 

292 

2S 

98 

,04 

.22 

.09 

-.18 

.80 

-12.6 

-18.7 

.05 

1106 

286 

75 

294 

.05 

3.02 

1.21 

-.51 

1.1 

•3.36 

-1.19 

-.17 

1111 

290 

75 

293 

.24 

1.98 

.80 

-.37 

.29 

-.74 

.55 

.22 

1121 

288 

climbing 

98 

,14 

-1.11 

-.44 

.90 

Summary 


•  The  combination  of  an  airborne  Doppler  Wind  Lidar 
and  a  Controlled  Towed  Platform  holds  promise  of  a 
transformation  of  how  we  Investigate  air/sea 
exchanges  and  construct  flux  parameterizations  for 
use  in  numerical  weather  models. 

•  This  research  is  being  funded  by  Dr.  Ferek  of  the 
Office  of  Naval  Research. 


Attachment  3 

Vertical  transport  of  TKE  within  Organized  Large  Eddies 
presented  at  the  DRI  workshop  2014 
by  G.  D.  Emmitt  (SWA),  S.  de  Wekker  (UVA)  and  R.  Foster  (Uwash) 


'  Vertical  transport  of  TKE  within 
Organized  large  Eddies 


G.  D.  Emmitt  (SWA),  S.  de  Wekker  (UVA)  and  R.  Foster  (Uwash) 

UPP  for  Extended  Range  Prediction 
DRI  Review 
Monterey,  CA 
27  -  28  August 


Objectives 


Investigate  and  characterize  the  presence  of 
rolls  (OLEs)  within  and  at  the  lateral 
boundaries  of  stratocumulus  topped  MBLs. 

Study  the  potential  impact  on  the 
development  and  implementation  of  the 
EDMF  into  forecast  models. 

Make  (or  not)  the  case  for  adding  an  MF  term 
unique  to  rolls  generated  by  shear  instabilities 
(non-convective). 


TODWL 

Twin  Otter  Doppler  Wind  Lidar 

CTV 

Controlled  Towed  Vehicle 


^SWA 


TODWL  data  description 


Attribute 

Performance  Metric 

Comments 

LOS  resolution  (applies 
to  vertical  profiles  of 

3D  winds  as  well) 

50  m 

Range  resolution  to  hard 
targets  (ground  or  dense 
cloud)  can  be  better  than 
10  meters. 

U,V,W  resolution 

<10  cm/s 

<  5  cm/s  for  stationary 
groundbased  operations 

Maximum  range 

6  -30  km 

Very  dependent  upon 
aerosols 

Time  to  complete  full 
step  stare  conical  scan 
for  wind  profiles 

-  20  sec 

12  point  step  stare  with 
.5  -2  second  dwells 

Sampling  frequency 

100  Hz 

Integration  of  several 
shots  is  typical  to 
improve  range 
performance 

Data  summary  for  September/October 
UPP  field  campaign  with  TODWL 

•  Three  quality  flights  with  the  DWL 

•  One  quality  co-flight  with  the  CTV  (9/30). 

•  Two  5  hour  ferry  flights  between  Monterey,  CA 
and  Salt  Lake  City,  Utah,  (Ruby  Mountain  waves) 

•  Seven  4.5  hour  flights  over  DPG,  Utah  during  the 
MATERHORN 

—  Co-funded  by  ONR  and  ARO 

—  3500  vertical  profiles  used  to  validate  WRF  models 

—  14  low  level  flights  dedicated  to  prospecting  for  OLEs 
in  the  vicinity  of  flux  towers. 


CIRPAS  Twin  Otter  with 
CTV  below 
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*  Derived  from  data  provided  by  Khelif  (UCI) 


structure  prospecting  with  DWL 


•  Feature  prospecting  uses  a  very  shallow  angle 
below  the  horizon  (~  -1  -3  degrees  for  a  300m 
flight  altitude). 

-  Results  in  ~  2  m  vertical  resolution  and  50  m 
horizontal  resolution  with  ~40  meter  sliding 
sample. 

-  It  takes  ~  33  seconds  to  profile  100  meters  below 
the  aircraft. 


TODWL  prospecting  geometries 

(9/30/12) 


300  m 


-3  degrees  elevation  angle  (tan  =  .052) 


E 


6  km 


Metric 

Value 

Twin  Otter  flight  altitude 

300  m  MSL 

Scan  configuration 

Straight  ahead  and  down  1  -3° 

Range  gate  length 

50  m 

Twin  Otter  ground  speed 

50  -60  m/s 

Shot  frequency 

160/second 

Shot  integration  (data  granule) 

100  (~  40  meters  forward  motion) 

Vertical  increment  between  granules 

~2  m 

Time  to  sample  100  m  vertical 

~  33  seconds 

Idealized  rolls 


Spectral  analyses  ofTODWLdata 


The  TODWL  is  a  2  micron,  eyesafe 
Coherent  detection  lidar.  The  returning 
Signal  is  heterodyned  with  a  local  optical 
Oscillator  and  then  digitized  (200  MHz). 

An  FFT  is  used  to  process  the  resulting 
Time  series  with  50m  line-of-sight 
Resolution.  In  homogeneous  flows,  a  single 
Spectral  feature  (top  figure)  is  identified 
as  the  frequency  representation  of  the 
radial  wind  speed. 

When  probing  a  turbulent  and/or  highly 
Sheared  environment,  the  spectral 
Feature  is  broadened.  It  is  this  broadening 
That  we  desire  to  use  for  estimating 
Turbulence  on  scales  of  25  -100  meters 
With  multipulse  integrations  over  50  -  200 
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9/30/12  case  study 

Processed  lidar  data  in  search  of  organized 
aerosol/wind  structures  below  the  Twin  Otter 
flight  path 

Processed  Twin  Otter  instrument  data  for  time 
series  of  u,  m,  w,  q,  and  0. 

Processed  CTV  data  for  u,v,w,q  and  0. 

Match  up  times  and  then  features  from  the 
TODWL  and  CABIN  data  sets  near  flight  level. 

Match  up  times  and  features  from  TODWL  and 
CTV  at  CTV  cruise  levels. 


Flight  path 


Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Mean  wind  at  fiight  levei  is  generaiiy  from  the  north 

Twin  Otter;  Racetrack  Twin  Otter:  Racetrack 
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Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Speed  (m/s) 

Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 
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Current  Findings 


Columns  of  ''spectral  broadening  (SB)"  seen 
only  in  zones  of  near  surface  convergence  and 
inferred  upward  motion 

Not  all  columns  terminate  at  elevated 
turbulence  "blobs"  or  accumulation  pockets. 

Number  of  SB  columns  per  1km  is  the  same  as 
the  number  of  convergence  zones.  However, 
the  spacing  between  convergence  zones 
varies. 


Line  density  of  rolls 
(per  10  km) 


Flight  Segment 

SB  Columns 
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Convergence  zones 
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CTV  Primes  102845  TKE:  0.80198  WQ: -12.655  W-Theta: -18.7389  Skew:  0.048147 
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Data  and  cross  products  derived  from  the  CTV  observations  (provided  by  Khelif  (CCTI) 
taken  during  the  case  above.  The  horizontal  axis  units  are  km  from  the  start  of  the  data  file. 
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Summary  of  segment  statistics 
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Summary  (Marine) 


The  combination  of  an  airborne  Doppler  Wind 
Lidar  and  a  Controlled  Towed  Platform  holds 
promise  of  a  transformation  of  how  we 
investigate  air/sea  exchanges  and  construct  flux 
parameterizations  for  use  in  numerical  weather 
models. 

Interpretation  of  the  TODWL  derived  spectral 
width  is  a  primary  focus  of  our  current 
investigation. 

Presentation  by  Foster  this  afternoon  will  address 
the  second  focus  area,  that  of  adding  a  "shear 
instability  driven  rolls"  term  to  the  EDMF. 


Rolls  over  land 


Extension  of  the  UPP  (Marine)  study  of  rolls 
contribution  to  the  EDMF  to  investigations 
over  land. 

Used  to  provide  comparisons  of  the  roll 
phenomena  as  measured  by  an  Airborne 
Doppler  Wind  Lidar. 

Flights  conducted  along  Salinas  Valley  and 
over  the  Dugway  Proving  Grounds,  Utah. 
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Box-and-whisker  plot  of  the  PBL  heights  derived  from  aerosol  backscatter 
profiles  along  the  north-south  flight  legs  (b)  performed  during  a  morning 
TODWL  mission  between  1012  and  1052  MDT  on  10  October  2012.  The 
horizontal  line  in  the  box  and  the  bottom  and  top  line  of  the  box  show  the 
median  of  the  data,  and  the  lower  and  upper  quartiles  (25%  and  75%), 
respectively.  The  whiskers  show  the  minimum  and  maximum  values  while 
the  solid  star  is  the  mean  value.  Vertical  potential  temperature  profiles 
obtained  over  the  Playa  and  Sagebrush  sites  are  shown  on  the  western  and 
eastern  part  of  the  domain.  Also  shown  is  the  potential  temperature  profiles 
obtained  from  DATAHAWK  (UAV)  measurements  located  between  TODWL 
legs  L2  and  L3. 
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Vertical  profiles  of  TODWL  derived  wind  speed  and  direction  along  a  north- 
south  flight  leg  east  of  Granite  Peak  (a)  and  an  east-west  flight  leg  over  Granite 
Peak  (b)  on  the  afternoon  of  14  October  2012. 
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From  top  to  bottom,  horizontal  profile  of  turbulence  kinetic  energy, 
vertical  velocity  vat  two  heights,  and  topography  below  a  flight  leg  east  of 
Granite  Peak  on  14  October  2012 


From  top  to  bottom:  w,  u,  and  v  -wind  component  along  a  flight  leg  east  of  Granite 
peak  on  14  October  2012.  The  flight  leg  shown  here  is  about  10  km  further  to  the 
east  than  the  flight  leg  shown  in  Fig.  3.  Shown  in  the  bottom  two  panels  are  the 
line  of  sight  velocity  and  the  signal-to-noise  ratio  of  a  forward  stare  on  the  same 
flight  leg.  In  these  panels,  coherent  turbulence  structures  are  visible  that  align  well 
with  velocity  fluctuations  from  the  in-situ  aircraft  measurements. 


Summary  (land) 


Both  convective  and  shear  driven  rolls  can  be 
seen  over  land  in  the  daytime. 

Salinas  Valley  exhibits  convective  rolls  in  the 
daytime  and  shear  driven  rolls  at  night  (not 
discussed  in  this  presentation). 

Deeper  boundary  layers  at  DPG  exhibit  rolls 
but  are  difficult  to  characterize  in  any  detail. 


Presentations  and  Students 


Two  graduate  students  supported  at  UVA. 
Focused  upon  the  DPG  data  sets. 

Papers  presented  at  multiple  working  groups, 
AGU  2013,  SPIE  2014,  AMS  2013,  2014. 

Paper  to  be  published  on  the  9/30  case  study 
on  hold  until  sampling  modeling  is  complete. 


Research  Plans 


Complete  DWL  sampling  model  developed  to  understand 
the  impact  of  various  scales  of  turbulence  and  shear  on 
spectral  broadening. (Emmitt  and  O'Handley) 

Reprocess  the  9/30/12  data  sets  guided  by  input  from  the 
sampling  model. (Emmitt  and  Godwin) 

Extend  Foster  model  to  cover  the  various  modes  of  the 
MBL  (Ri,  shear,  convective  stability)  (Foster) 

Conduct  sensitivity  experiments  with  a  modified  EDMF 
within  the  WRF  and  perhaps,  the  COAMPS,  (de  Wekker  and 
Foster). 

Produce  a  lessons  learned  document  and  draft  plan  for  an 
intensive  TODWL/CTV  investigation  of  non-convective  rolls. 
Interact  with  the  roll  investigators  (Mahrt,  Ginis, ..)  in  the 
design. 
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*  Derived  from  data  provided  by  Khelif  (UCI) 


Attachment  4 

EDMF  in  the  Near-Neutrally  Stratified  Boundary  Layer 
presented  at  the  DRI  workshop  2014 
by  R.  Foster  and  G.  D.  Emmitt 


EDMF  in  the  Near-Neutrally 
Stratified  Boundary  Layer 


Ralph  Foster,  APL  UW 
G.  David  Em  mitt,  SWA 
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Cold  air  flow  over  ice: 
smaller-scale  rolls, 
near-neutral,  possibly 
slightly  stable 


Cold  air  flow  off  ice, 
over  warmer  water: 
larger-scale  rolls,  merging. 
Broadening,  deepening 


JC.  -S. 

iif'*  ^ 


^  S'  -  ,#1'' 

PBLtransition  ^ nn” 
roll  regirnelnto;;  * 
CO n ve ct i ve  regi  m e  > 

<  “  j;-, 


V 


EDMF 


Method  for  estimating  the  non-gradient 
contribution  to  the  PBL  fluxes 

Current  MF  focus  is  on  non-gradient  contribution 
in  low  winds,  primarily  buoyancy  production 

-  Assumes  ED  model  correctly  models  PELS  with 
significant  shear  production 

Does  EDMF  need  to  capture  the  non-gradient 
fluxes  due  to  roll  OLE? 

—  Extension  of  MF  contribution  to  near-neutral 
stratification 


PBL  Rolls 


Very  common 

-  Span  the  depth  of  the  PBL 

-  Due  to  combined  shear  &  convective  instability 

-  UlO  >~  5  to  7  m/s 

-  Near-neutral  to  moderately  unstable  stratification 

-  Non-gradient  diffusion  transport 

•  Is  this  contribution  important? 

•  Can  we  measure  it? 

•  Can  we  parameterize  it? 
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ill  the  hurricane  boundary  layer  by  large  eddies 

Ping  Zhu‘ 

Received  27  November  2007:  revised  "'I  Mardi  2008;  accepted  23  May  2008;  published  3  September  2t<)8 

(i]  Hurricane  boundaiv  layer  (I113L)  proces.scs,  especially  the  stnictuie  of  the  coherent 
lat^c  eddy  circulatioii.s  (LI  .C.s)  and  tiieir  induced  vertical  tran.sport,  are  not  well 
understood.  This  paper  intnxluccs  a  large  eddy  siimilation  (LliS)  t’rantework  in  a  weatliei 
hindcasting  mode  developed  trom  a  nmltiplc  scale  nested  Weather  Research  and 
l  orecasiing  (WRF)  model.  Using  the  WRF'-LLS.  this  study  investigated  the  slmclurc  of 
the  IIBL  LL(.*s  and  the  a.ssociaied  vertical  transpon  during  the  landfall  of  Hurricane  I\an 
(2004).  The  simulation  slums  that  the  IIBL  LliCs  exist  in  a  mean  stable  envimiunent  and 
consist  of  vvell'dcrmcd  updraft  and  downdralT  Statistically,  the  llBL  LLC's  are  only 
sligjilly  skewed  with  (lie  iipdratls  and  downdmlls  ivlali\cly  evetily  distributed  spatially. 


(c)  Large  eddy  circulations  and  virtual  potential  temperature  (k)  at  the  cross  section  in  (a.3) 
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(d)  Large  eddy  circulatiwis  and  mixing  ratio  (g  kg“’)  at  the  cross  section  in  (a.3) 


X~distanco  (km) 


4  ZHU;  HURRICANE  BOUNDARY  LAYER  TURBULENT  TRANSPORT 

Distance  between  the  cores  of  maximum  and  minimum  winds  (km) 
30 1 - ^ ^ - 1 - 1 - 1 - 1 - 1 - 1 - r 


25 


Figure  5.  Histograin  of  the  di.stance  between  the  cores  of  maximum  and  minimum  winds  in  the  surface 
wind  bands. 
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Figinv  fi,  Flux  profiles  of  resolved  turbulent  eddies  and  SOS  turbulence  averaged  over  the  laidfai) 
period  fium  05:(X)  UTC  to  08:00  UTC.  Bladj  solid:  total  resolved  fluxes;  daslied:  fluxes  induced  by  the 
IIBL  LECs;  dotted:  fluxes  generated  by  chaotic  hffbulcnt  eddies.  Light  solid:  fluxes  of  S<iS  turbulence. 
The  height  has  been  normalized  by  the  mean  HBL  height,  which  is  defined  as  the  base  of  the  inversion 
layer 
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Figure  12.  (a)  Skewness  of  LB(.'s  averaged  over  the 
landfall  period  from  05:00  UTC  to  08.00  UTC.  (h)  The 
mean  spatial  fraction  of  updraft  over  the  same  period. 


I'igurv  9,  The  resolved  and  S(iS  Irrbulent  fluxes  in  D-5  avenged  over  the  land&Il  period  fiom 
05:00  UTC  m  08:(K)  UTC  compared  with  the  parameterized  turbulent  fluxes  at  Ihe  besrtion  in  D-L  0-2, 
and  D3  where  D-5  is  nested.  The  dark  solid  and  Ihe  light  dot-dashed  profiles  indicate  the  total  resolved 
and  S(iS  fluxes  in  D-5,  respectiv'ely.  The  light  thick  dashed  profile  indicates  the  contributirwi  from  the 
resolved  chaotic  eddies  only.  The  light  solid,  light  dotted,  and  dark  dasfied  profiles  arc  the  parameterized 
turbulent  fluxes  from  domains  D-l,  D-2,  and  D-5  at  the  grid  pwint  where  D-5  is  nested,  respectively. 


5-domain  WRF  "LES"  (D5:  ~100  x  10  m) 


Last  UPP  Meeting 


It  appears  that  a  PBL  jet  can  induce 
"stacked"  ole 

Co-linear  OLE  triads  can  generate 

COMPLEX  STRUCTURE 

There  is  a  weaker,  larger-scale,  oblique 

ABOVE-JET  MODE 

Can  it  be  energized  by  nonlinear 
triad  mechanism? 

Future  work 

Modify  triad  code  for  non-co-linear 
OLE 

(2014:  Done,  preliminary  look 
shows  highly  complex  roll  OLE) 
Include  stratification  and  thermal 
wind  effects  on  OLE 

(2014:  Done,  untested) 


lH"  (shading)  anU  i|r 


=  5.03  (ND):  Xp  =  7.85  (ND);  =13.96  (ND) 

W  (shading)  and  y  (contour):  time  =  0150  (ND) 
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Rolls  in  EDMF  Framework 


Simple  MBL  conditions 
Look  at  basic  co-linear  triads 
Simple  MF  statistics 

-  w'  Skewness 

-  Updraft  fraction  [a] 

-  Large-eddy  mass  flux  {a(l  -  o-)(w„p  -  w^n)) 

Compare  with  previous  work  (Zhu,  2007) 
Compare  with  Observations 


Mean  Flow  U(black)  V(red)(m/s) 


Basic  PBLflow 
(Compared  to 
pure  Ekman) 


0.06 

0.05 

0.04 

0.03 

0.02 

0.01 

0 

^-0.01 

-0.02 

-0.03 

-0.04 

-0.05 

-0.06 


=  O.b'l  =  '\Z.O/  {WUy,  =  lv3.yD  {\\U) 

W  (shading)  and  vj;  (contour):  time  =  1275  (ND) 


U"^  (shading)  and  \\f  (contour) 


5  10  15  20 

x(ND) 


Z(ND) 


Updraft  Fraction  (a)  w'  Skewness 


0.6 


For  this  triad,  statistics  vary  little  with  time 


w'v 


/ 


Discussion 


Idealized  (nonlinear,  triad  interaction)  roll  OLE 
model  produces  results  qualitatively  similar  to  WRF- 
"LES"  of  hurricane  PBL  LEC  effects 

-  Similar  LEC  skewness  and  updraft  fraction 

—  LEC  contribution  much  larger  than  PBL  parameterization 

Differences  with  Monterey  Bay  lidar  observations 

—  Flux  data  suggested  negative  skewness  near  top  of  OLE 

-  Obs  suggest  OLE  effect  on  fluxes  not  as  large  as  modeled 

•  (10  to  30%,  not  half  or  more) 

•  Lateral  mixing  missing  in  model? 

•  Direct  OLE  effect  may  overestimate  contribution;  "collection"  of 
turbulent  eddies  may  be  more  important 

•  Low-order  truncation  exaggerates  fluxes 


Low-Order  Truncation  Errors 


Re  =  500.  alpha  =  0.500.  epsilon  —  16.00  deg 


Mean-flow  Modification 


Re  =  150.  alpha  =  0.540.  epsilou  =  1 1 .00  deg 


Amplitude  Estimation 


0 

0 

0 

0, 

0 

0, 

0 


Low-order  truncation  problems: 

•  Over-estimated  amplitude 

•  "S-shaped"  MF  modification 
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Standard  Non-Linear  Single-Wave  PBL  Roll  Model 


Table  5.1  Coutribiitioiis  to  the  iioiiUuear  perturbation  up  to  the  fifth  Lauclau 
Coefficient 


Truncated  Contributions 
to  Multi-Wave  Roll  Model 


q=  [u,  V,  w,  T]t 


0  + 

Aq^o  +  0  + 

0  +  A^q2o 

0  +0 


(mean  flow  modification) 
(fundamental  wavelength) 
(1^^  harmonic) 

(2"^  harmonic) 


+  A^qjo 
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Triad  Interaction  Model 


qa  = 


qp 


qv  = 


+  c[A^q2.v  ^^q3,v  ^  ^^q4.v] 


14, V- 


:  contributions  from  single-wave  theory;  e.g.  q2  „  =  „  • 

•  Blue:  new  wave-wave  &  wave-mean  flow  interaction  contributions. 

•  Red:  Low-order  phase-coupling  terms. 

•  Also:  mean-flow  modifications  due  to  each  wave. 


15 


Z  (ND)  Z  (ND) 


W  (shading)  and  v(/  (contour);  time  =  1250  (ND) 


W  (shading)  and  \|/  (contour);  time  =  1350  (ND) 


(shading)  and  \\f  (contour) 


1-0.1 
-02 
-O.’S 


U'^  (shading)  and  \\f  (contour) 


Similar  triad  generates  different  OLE  effect 


V 

-8.6  -0.4 


w’  Skewness 


More  similar  to  Obs, 

Less  similar  to  WRF-"LES'' 


For  this  triad,  statistics  do  vary  with  time 
Hurricane  is  very  high  shear,  ideal  for  roll  OLE, 
strongly  forced  modes. 

Monterey  Bay  had  more  variable  conditions  & 
layered  PBL  structure. 


^  ^  &  e  P 


What  Do  We  Need  From 
Observations? 


Need  to  constrain  parameter  space 

—  PBL  profiles  throughout  roll  prospecting  region 

—  For  given  mean  flow  (U,VXq),  what  are  the 
characteristics  of  the  rolls?  {wavelength(s), 
orientation,  strength) 

—  Can  we  characterize  in  terms  of  Ri,  Re, ...? 

Effects  of  rolls  on  fluxes: 

—  Fluxes  and  flux  statistics  at  multiple  levels 

•  Sfc,  mid-PBL,  Upper  levels 

•  Separate  local  from  LECflux  contributions 


Next  Steps 


Investigate  stratification  effects  &  3-D  effects 

Use  observations  to  correct  model 

Implement  simple  roll  MF  addition  to  EDMF 
model 


Attachment  5 

Investigation  of  coastal  zone  complex  boundary  circulations  and  their  potential  impact  on  energy 

transport 

presented  at  SPIE  Defense  and  Security  Conference  in  Amsterdam  2014 
by  G.  D.  Emmitt,  K.  Godwin  and  R.  Foster. 


Investigation  of  coastal  zone 
complex  boundary  circulations  and 
their  potential  impact  on  energy 

transport 


G.  D.  Emmitt  and  K.  Godwin 
Simpson  Weather  Associates,  Inc. 

R.  Foster 

University  of  Washington 


SPIE  Amsterdam  2014 


Objectives 


Investigate  and  characterize  the  presence  of  rolls 
(OLEs)  within  and  at  the  lateral  boundaries  of 
stratocumulus  topped  Marine  Boundary  Layers 
(MBLs). 

Study  the  potential  impact  on  the  development 
and  implementation  of  the  Eddy  Diffusivity  and 
mass  Flux  (EDMF)  parameterization  into  weather 
forecast  models. 

Make  the  case  for  adding  an  MF  term  unique  to 
rolls  generated  by  both  convective  and  shear 
instabilities. 


CTV 


Surface 
Temperature 
Sensor 


TODWL 

Twin  Otter  Doppler  Wind  Lidar 

CTV 

Controlled  Towed  Vehicle 


TODWL 

scanner 


SWA 


TODWL  data  description 


Attribute 

Performance  Metric 

Comments 

LOS  resolution  (applies 
to  vertieal  profiles  of 

3D  winds  as  well) 

50  m 

Range  resolution  to  hard 
targets  (ground  or  dense 
cloud)  can  be  better  than 
10  meters. 

U,V,W  resolution 

<10  cm/s 

<  5  cm/s  for  stationary 
groundbased  operations 

Maximum  range 

6  -30  km 

Very  dependent  upon 
aerosols 

Time  to  eomplete  full 
step  stare  eonieal  sean 
for  wind  profiles 

~  20  sec 

12  point  step  stare  with 
.5  -2  second  dwells 

Sampling  frequeney 

100  Hz 

Integration  of  several 
shots  is  typical  to 
improve  range 
performance 

Data  summary  for  September/October 
DPP  field  campaign  with  TODWL 

•  Three  quality  flights  with  the  DWL 

•  One  quality  co-flight  with  the  CTV  (9/30). 

•  Two  5  hour  ferry  flights  between  Monterey,  CA 
and  Salt  Lake  City,  Utah.  (Ruby  Mountain  waves) 

•  Seven  4.5  hour  flights  over  DPG,  Utah  during  the 
MATERHORN 

-  Co-funded  by  ONR  and  ARO 

-  3500  vertical  profiles  used  to  validate  WRF  models 

-  14  low  level  flights  dedicated  to  prospecting  for  OLEs 
in  the  vicinity  of  flux  towers. 


CIRPAS  Twin  Otter  with 
CTV  below 


Cabin  Primes  102845  TKE:  0.038198  WQ:  0.21743  W-Theta:  0.087407  Skew; -0.1774 

U-prim*  V*p(im«  W-pnm* 


TTvxBHn  bvH 


Twin  Otter 
Probe  Data 


CTV  Primtt  102845  TKE:  880188  Wa>1L#5S  -18.758a  £kew:  fl  04B147 
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Flux  Data* 


*  Derived  from  data  provided  by  Khelif  (UCI) 


structure  prospecting  with  DWL 


•  Feature  prospecting  uses  a  very  shallow  angle 
below  the  horizon  (~  -1  -3  degrees  for  a  300m 
flight  altitude). 

-  Results  in  ~  2  m  vertical  resolution  and  50  m 
horizontal  resolution  with  ~40  meter  sliding 
sample. 

—  It  takes  ~  33  seconds  to  profile  100  meters  below 
the  aircraft. 


TODWL  prospecting  geometries 

(9/30/12) 


-3  degrees  elevation  angle  (tan  =  .052) 


6  km 


Metric 

Value 

Twin  Otter  flight  altitude 

300  m  MSL 

Scan  configuration 

Straight  ahead  and  down  1  -3° 

Range  gate  length 

50  m 

Twin  Otter  ground  speed 

50  -60  m/s 

Shot  frequency 

160/second 

Shot  integration  (data  granule) 

100  (~  40  meters  forward  motion) 

Vertical  increment  between  granules 

~2  m 

Time  to  sample  100  m  vertical 

~  33  seconds 

Idealized  rolls 


Convergence 


Divergence 


Spectral  analyses  of  TODWL  data 


The  TODWL  is  a  2  micron,  eyesafe 
Coherent  detection  lidar.  The  returning 
Signal  is  heterodyned  with  a  local  optical 
Oscillator  and  then  digitized  (200  MHz). 

An  FFT  is  used  to  process  the  resulting 
Time  series  with  50m  line-of-sight 
Resolution.  In  homogeneous  flows,  a  single 
Spectral  feature  (top  figure)  is  identified 
as  the  frequency  representation  of  the 
radial  wind  speed. 

When  probing  a  turbulent  and/or  highly 
Sheared  environment,  the  spectral 
Feature  is  broadened.  It  is  this  broadening 
That  we  desire  to  use  for  estimating 
Turbulence  on  scales  of  25  -100  meters 
With  multipulse  integrations  over  50  -  200 


C\UPP  2012>.DASi.093020i2i_100709s0^ 


CMJPP_  20i2\DAS\0930201 2{_^0070960i 


9/30/12  case  study 

Processed  lidar  data  in  search  of  organized 
aerosol/wind  structures  below  the  Twin  Otter 
flight  path 

Processed  Twin  Otter  instrument  data  for  time 
series  of  u,  m,  \n,  q,  and  0. 

Processed  CTV  data  for  u,v,w,q  and  0. 

Match  up  times  and  then  features  from  the 
TODWL  and  CABIN  data  sets  near  flight  level. 

Match  up  times  and  features  from  TODWL  and 
CTV  at  CTV  cruise  levels. 


Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


17:00 


Z{m) 


Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Flight  path 


Latitude 


0930  Flight  Path  -  1 00709  Mean  Heading;  93 


LOSVelocity  093012  100709  Az.  Range: -8  to  2  El.  Range: -3  to  0 
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Spectral  Width 
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0930  Flight  Path  -  101709  Maan  Heading:  299 
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0930  Flight  Path  - 1 1 0626  Mean  Heading:  294 
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Spectral  Width 
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093012-102845  Heading:  98.1174 


LOS  Velocity  093012  102845  Az.  Range: -7  to -6  El.  Range: -3  to -2 


Distance  (km) 


Current  Findings 


Columns  of  "spectral  broadening  (SB)"  of  the 
DWL  return  are  seen  only  in  zones  of  near  surface 
convergence  and  inferred  upward  motion 

Not  aii  columns  terminate  at  elevated  turbulence 
"blobs"  or  accumulation  pockets. 

Number  of  SB  coiumns  per  1km  is  the  same  as 
the  number  of  convergence  zones.  However,  the 
spacing  between  convergence  zones  varies 
greatiy. 


Line  density  of  rolls 
(per  10  km) 


Flight  Segment 
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Summary 


The  combination  of  an  Airborne  Doppler  Wind 
Lidar  and  a  Controlled  Towed  Platform  holds 
promise  of  a  transformation  of  how  we 
investigate  air/sea  exchanges  and  construct  flux 
parameterizations  for  use  in  numerical  weather 
models. 

Interpretation  of  the  TODWL  derived  spectral 
width  is  a  primary  focus  of  our  current 
investigation. 

We  are  now  focused  on  adding  a 
"convective/shear  instability  driven  rolls"  terms 
to  the  EDMF. 


Attachment  6 


Triad  Interaction  Model  and  its  role  in  modifying  the  EDMF  parameterization  to  include 
contributions  from  non-gradient  fluxes  due  to  rolls. 

Sub  report 

R.  Foster 


Triad  Interaction  Model  and  its  role  in  modifying  the  EDMF  parameterization  to  include 
contributions  from  non-gradient  fluxes  due  to  rolls. 


Ralph  Foster 


An  important  part  of  our  project  is  to  use  the  Doppler  wind  lidar  (DWL)  data  to  guide  in  the  development  of  a 
parameterization  of  roll  OLE  contributions  to  the  PBL  fluxes.  The  target  for  this  parameterization  is  the  eddy 
diffusivity/Mass  flux  (EDMF)  PBL  parameterization.  As  currently  implemented,  the  EDMF  parameterizes  the 
non-gradient  fluxes  that  are  induced  by  convective  plumes  using  a  mass  flux  scheme.  The  plume  contribution 
parameterization  is  based  LES  of  shear-free  convection.  However,  when  shear  is  important,  and  even  into  the 
very  near-neutral  stratification  regime,  the  dominant  form  of  non-gradient  fluxes  is  due  to  roll  OLE,  which  are 
much  more  challenging  to  study  and  parameterize  using  LES. 

For  example,  an  EDMF  model  must  be  capable  of  capturing  complex  PBL  transitions  such  as  in  off-ice  flow  in 
high  latitudes.  The  PBL  is  slightly  stable  over  pack  ice,  yet  rolls  are  observed  in  this  regime.  Over  water,  that 
flow  gains  significant  surface  buoyancy  flux  and  the  rolls  strengthen,  broaden  and  deepen.  Farther  off-shore, 
they  experience  an  up-scale  transition  and  ultimately  break  down  into  0(25-50  km)  convection.  At  present,  the 
EDMF  parameterizes  the  small-scale,  largely  shear-driven  small-scale  turbulence  with  the  ED  part  and  the 
convective  regime  after  the  rolls  break  down.  However,  in  the  near  neutral  stable  and  the  moderately  unstable 
regimes  where  rolls  are  present,  the  EDMF  model  has  no  mechanism  for  parameterizing  their  non-gradient 
contribution  to  the  PBL  fluxes. 

Our  strategy  is  to  use  the  DWL  observations  to  calibrate  and  validate  a  general  model  for  roll  OLE  that  can  be 
used  to  develop  a  roll  parameterization  for  the  EDMF.  Roll  OLE,  which  are  often  mischaracterized  as  either 
convective  rolls  or  shear  rolls  are  in  general  due  to  a  combined  shear  and  convective  instability  that  reaches  a 
quasi-equilibrium  state  in  which  finite  amplitude  rolls  are  maintained  in  a  slightly  modified  mean  state.  This 
mechanism  is  well-established  in  the  literature,  most  recently  for  the  hurricane  boundary  layer  by  Foster 
(2005),  Zhu  (2008)  and  Nakanishi  and  Niino  (2012).  However,  both  our  DWL  observations  and  the  numerical 
modeling  indicate  that  single-mode  rolls,  while  capturing  the  basic  OLE  flow,  are  insufficient  to  describe  the 
variability  in  the  PBL. 

As  previously  reported,  the  general,  nonlinear  "single-mode"  model  described  in  Foster  (2005)  has 
been  expanded  to  capture  the  nonlinear  energy  transfer  among  resonant  triads  of  PBL  modes.  Resonant  triads 
are  the  fundamental  building  blocks  of  the  nonlinear  interactions  between  fluid  motions  at  different  scales. 
Energy  transferred  from  motions  of  different  spatial  scales  into  or  out  of  a  particular  mode  comes  from  modes 
whose  wave  vectors  form  a  triangle  with  the  wave  vector  for  the  mode  in  question  (Fig.  1).  The  standard 
single-mode  nonlinear  stability  analysis  (e.g.  Foster.  2005)  will  thus  omit  these  "wave-wave  interaction" 
pathways  by  which  energy  can  flow  from  the  dominant  mode  into  more  slowly-growing  modes.  Our  first  effort 
worked  on  the  much  simpler  case  in  which  the  triad  modes  are  all  colinear. 


In  this  year,  we  significantly  expanded  the  model  to  include  the  nonlinear  interactions  between  non-colinear 
triads.  The  quasi-equilibrium  solutions  continue  to  change  with  time,  but  through  a  limited  flow  regime.  And  at 
any  particular  time,  the  solutions  vary  spatially.  The  solution  space  is  complex,  but  so  far  the  results  are 
realistic  and  appear  to  agree  with  simpler  models  (e.g.  Dubos  et  al.,  2008).  A  snapshot  of  a  particular  solution 
is  shown  in  Figs.  2  and  3.  The  dominant  roll  OLE  are  oriented  about  19°  from  the  down-wind  direction.  The 
triad  includes  a  smaller-scale  mode  oriented  about  17°  from  the  dominant  mode  and  a  larger-scale  mode 
about  46°  from  the  dominant  mode.  The  overturning  circulation  induced  by  the  roll  OLE  is  best  seen  in  the  x-z 
cross-sections  taken  at  different  down-roll  positions  (Fig.  2).  The  variability  in  the  overturning  circulation 
strength  and  shape  is  evident,  which  can  be  seen  most  clearly  in  the  vertical  velocity.  The  along-roll  velocity 
also  varies  in  strength  along  the  rolls.  Fig.  3  shows  the  horizontal  perturbation  wind  vectors  at  near  the  surface 
and  in  the  mid-PBL.  near  the  surface,  the  along-roll  component  dominates.  In  the  mid-PBL,  the  overturning 
circulation  dominates.  The  roll  circulation  is  clearly  three-dimensional  and  varies  in  both  time  and  space. 
Efforts  in  this  year  will  focus  on  understanding  the  non-colinear  triad  interactions  in  this  simple  case.  This 
study  is  an  important  advance  in  PBL  roll  theory  and  is  worthy  of  publication.  Once  the  new  (and  more 
complex)  mathematics  is  more  carefully  checked  and  better  understood,  it  will  be  applied  to  the  DWL 
observations. 

Zhu  (2008)  proposed  that  roll  OLE  fluxes  could  be  parameterized  using  the  mass-flux  paradigm.  A  key 
difference  is  that  standard  mass  flux  modeling  assumes  that  the  vertical  velocity  is  highly  skewed  with  a  small 
updraft  fraction.  For  roll  OLE,  the  vertical  velocity  is  only  slightly  skewed  and  the  updraft  fraction  is  closer  to 
50%.  For  a  single  roll  OLE  situation,  the  skewness,  updraft  fraction  and  mass  flux,  (l-a)(w_up-w_dn  ),  are  all 
well  behaved  and  similar  to  the  calculations  in  Zhu  (2008)  (Fig.  4  and  5).  As  shown,  OLE  roll  fluxes  are 
inconsistent  with  eddy-diffusivity  modeling.  However,  the  mass  flux  parameterization  is  relatively  well- 
behaved. 

In  support  of  the  application  to  the  EDMF,  we  applied  the  colinear  triad  interaction  model  to  the  OLE 
observations  over  Monterey  Bay.  Previously,  we  made  a  case  that  the  mean  PBL  structure  was  likely  due  to 
advection  of  over-land  PBL  air  over  the  thin  marine  layer  in  the  bay.  The  mean  flow  indicated  an  elevated  PBL 
jet  at  the  intersection  between  the  air  masses.  The  observations  showed  a  likely  "stacked-OLE"  configuration 
in  which  OLE  were  separately  generated  in  each  layer  and  interacted  across  the  jet.  The  triad  model  showed 
greater  and  lesser  interaction  between  the  layers  above  and  below  the  jet  that  varied  in  spatially  and 
temporally,  which  was  consistent  with  the  observations  (Fig.  6).  Note  the  change  in  the  Skewness  profiles 
above  the  level  of  the  mid-PBL  jet. 

In  this  case,  consistent  with  the  DWL  data.  Twin  Otter  fluxes  and  towed  platform  fluxes,  the  OLE-induced 
skewness  is  much  more  complex.  Understanding  and  modeling  correctly  this  process  will  be  a  focus  in  the 
upcoming  year.  In  particular,  we  will  need  to  implement  the  non-colinear  roll  OLE  model  and  apply  it  to  the 
Monterey  PBL  conditions. 

Joao  Teixeira's  research  group  at  JPL  has  provided  us  with  an  implementation  of  the  EDMF  PBL  model  that  we 
can  use  to  experiment  with  implementing  the  OLE  MF  contribution. 


Figure  1:  Typical  non-colinear  roll  OLE  triad  for  simple  PBL  mean  flow.  Contours  show  the  instability  growth 
rate.  Dominant  mode  is  labeled  B.  Sub-modes  are  A  and  C.  Construction  shows  that  A  =  B  +  C  as  required. 

Figure  2;  Cross-sections  in  the  x-z  plane  at  two  locations  along  the  dominant  roll  mode  direction  (mode  B  from 
Fig.  1).  Upper  panels  color  shading  show  vertical  velocity.  Lower  panels  color  shading  shows  along-mode  B  OLE 
perturbation  velocity.  Vectors  show  the  overturning  circulation  (U,  W)  perpendicular  to  the  mode  B  OLE  axes. 

Figure  3:  Horizontal  cross  sections  (x-y)  through  the  vertical  velocity.  Left  panel  is  nearer  the  surface  and  right 
panel  is  at  the  mid-PBL.  Vectors  show  the  along-across  horizontal  velocities.  Near  the  surface,  the  along-mode 
B  component  dominates.  In  the  mid-PBL,  the  overturning  component  begins  to  dominate.  The  secondary 
circulation  due  to  rolls  is  three-dimensional. 

Figure  4:  Upper  panel  pair:  Snapshot  of  colinear  roll  OLE  triad  solution  for  a  simple  PBL  mean  flow.  Upper 
shows  vertical  velocity  and  overturning  stream  function.  Lower  shows  along-roll  perturbation  velocity.  Lower 
left  panel:  updraft  fraction.  Lower  right  panel:  vertical  velocity  skewness. 

Figure  5:  Example  across-roll  momentum  flux.  Upper  panel,  wu  velocity  product  with  overturning 
streamfunction  contoured.  Lower  left:  effective  eddy  viscosity  needed  if  these  fluxes  are  modeled  using  the 
eddy  diffusivity  approach.  Lower  right,  momentum  flux  modeled  using  mass  flux  approach. 

Figure  6:  Upper  panels:  two  snap  shots  of  the  vertical  velocity  for  the  nonlinear,  colinear  triad  model  for  a 
mean  flow  profile  similar  to  the  DWL  observations  over  Monterey  Bay,  2012.  Lower  left:  mean  flow  profile 
used  to  generate  the  solution  show  in  the  upper  panels.  Lower  right:  vertical  velocity  skewness  from  the  triad 
model.  Note  that  the  skewness  can  be  either  positive  or  negative  near  the  surface,  which  is  consistent  with  the 
observations. 
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Investigation  of  the  representation  of  OLEs  and  terrain  effects  within  the  coastal  zone  in  the 
EDMF  parameterization  scheme:  an  airborne  Doppler  wind  lidar  perspective. 

Stephan  FJ.  De  Wekker 
University  of  Virginia 

Reporting  period:  1  July  2013-30  June  2014 

The  planned  tasks  for  year  3  of  this  project  that  involved  subcontractor  UVA  were: 


Task  3.1  -  Processing  and  analysis  (including  algorithm  and  software  development) 
of  new  TODWL  data  after  completion  of  field  program 

Task  3.3  -  Analysis  of  turbulence  obtained  from  the  standard  instruments  on  the 
Twin  Otter 

Task  3.7  -  Analyze  TODWL  data  and  compare  with  Radiosonde  data  (UVA  and  SWA) 
Accomplishments: 

During  year  3,  a  majority  of  the  effort  was  in  the  analysis  of  TODWL  research  flights.  The  UVA 
contribution  was  focused  on  investigating  of  turbulence  and  boundary  layer  structure  over  land,  in 
particular  during  the  ONR  funded  MATERHORN  field  campaign.  Seven  research  flights  in  total  were 
conducted  during  the  MATERHORN  field  campaign. 

We  analyzed  PEL  heights  from  the  TODWL  data  for  all  research  flights.  An  example  of  the  results  for 
five  flight  legs  oriented  north-south  for  the  morning  of  10  October  2012  is  shown  in  Fig.  1.  A  terrain 
following  PEL  height  is  evident  with  maximum  heights  up  to  about  2000  m.  The  TODWL  derived 
PEL  heights  compare  well  with  PEL  height  derived  from  temperature  profiles  from  radiosonde  and 
unmanned  aerial  vehicle  measurements. 

A  typical  feature  of  the  boundary  layer  structure  around  Granite  Peak  is  the  presence  of  layers  with 
large  amounts  of  shear  in  wind  speed  and  direction.  We  hypothesize  the  presence  of  these  layers  play 
an  important  role  in  modulating  turbulence  intensities  and  turbulence  structures.  An  example  is  shown 
for  the  afternoon  of  14  October  2012  for  a  north-south  oriented  flight  leg  east  of  Granite  Peak.  In  a 
layer  below  the  PEL  height  at  around  2200  m,  there  was  large  shear  in  wind  speed  and  direction, 
especially  on  the  northeastern  part  of  Granite  Peak  (Fig.  2).  In  situ  measurements  from  the  Twin  Otter 
aircraft  indicate  enhancement  of  turbulence  kinetic  energy  on  the  northern  part  of  this  flight  leg  (Fig. 
3).  The  low  level  Twin  Otter  flights  were  designed  to  prospect  for  organized  convective  structures  and 
for  quantifying  vertical  heat  and  momentum  fluxes  associated  with  these  structures.  Forward  stares 
from  the  lidar  indicate  the  presence  of  organized  structures  that  align  well  with  the  in-situ  turbulence 
measurements  (Fig.  4) 
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Figure  1:  (a)  Box-and-whisker  plot  of  the  PBL  heights  derived  from  aerosol  backscatter  profiles  along  the  north-south  flight  legs  (b) 
performed  during  a  morning  TODWL  mission  between  1012  and  1052  MDT  on  10  October  2012.  The  horizontal  line  in  the  box  and  the 
bottom  and  top  line  of  the  box  show  the  median  of  the  data,  and  the  lower  and  upper  quartiles  (25%  and  75%),  respectively.  The  whiskers 
show  the  minimum  and  maximum  values  while  the  solid  star  is  the  mean  value.  Vertical  potential  temperature  profiles  obtained  over  the 
Playa  and  Sagebrush  sites  are  shown  on  the  western  and  eastern  part  of  the  domain.  Also  shovm  is  the  potential  temperature  profiles 
obtained  from  DATAHAWK  (UAV)  measurements  located  between  TODWL  legs  L2  and  L3. 
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Figure  2:  Vertical  profiles  of  TODWL  derived  wind  speed  and  direction  along  a  north-south  flight  leg  east  of  Granite  Peak  (a)  and  an 
east- west  flight  leg  over  Granite  Peak  (b)  on  the  afternoon  of  14  October  2012. 
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Figure  3:  From  top  to  bottom,  horizontal  profile  of  turbulence 
kinetic  energy,  vertical  velocity  vat  two  heights,  and 
topography  below  a  flight  leg  east  of  Granite  Peak  on  14 
October  2012 


Figure  4:  From  top  to  bottom:  w,  u,  and  v  -wind  component 
along  a  flight  leg  east  of  Granite  peak  on  14  October  2012.  The 
flight  leg  shown  here  is  about  10  km  further  to  the  east  than  the 
flight  leg  shown  in  Fig.  3.  Shown  in  the  bottom  two  panels  are 
the  line  of  sight  velocity  and  the  signal-to-noise  ratio  of  a 
forward  stare  on  the  same  flight  leg.  In  these  panels,  coherent 
turbulence  structures  are  visible  that  align  well  with  velocity 
fluctuations  from  the  in-situ  aircraft  measurements. 


These  coherent  turbulence  structures  are  sub-grid  scale  for  regional  to  global  scale  models  and 
therefore  one  would  expect  an  EDMF  based  parameterization  scheme  to  represent  the  effects  of  these 
structures  on  the  turbulent  transport.  We  are  currently  collaborating  with  Dr.  James  Doyle  at  NRL  in 
the  use  of  the  CO  AMPS  to  simulate  these  MATERHORN  case  studies.  Doyle’s  group  has  now 
implemented  the  EDMF  scheme  in  COAMPS  and  will  perform  simulations  for  one  of  the  intensive 
operational  periods  during  MATERHORN.  Simulations  will  be  done  with  and  without  EDMF 
parameterization  at  horizontal  grid  spacing  of  4,3,2,  and  1  km  in  the  innermost  domain  with  the  overall 
goal  to  evaluate  the  sensitivity  of  the  simulations  to  ED  and  MF  parameters  in  the  EDMF 
parameterizations.  PBL  structure  and  the  simulated  fluxes  of  heat  and  momentum  in  the  PEL  will  be 
compared  against  the  in-situ  and  lidar  measurements  from  the  Twin  Otter  flights.  In  addition,  De 
Welder  has  continued  a  collaboration  with  the  Research  Application  Laboratory  at  the  National  Center 
for  Atmospheric  Research  (NCAR-RAL  -  Knievel  and  Liu)  to  perform  ‘very  large  eddy  simulations’ 
(VLES)  with  a  300  m  grid  spacing  for  the  various  days  of  the  research  flights.  Idealized  simulation 
have  also  been  performed  in  collaboration  with  the  University  of  Vienna  (Serafln).  Output  will  in  the 
next  year  be  analyzed  to  detect  the  presence  of  coherent  structures  similar  as  in  the  observations  and 
quantify  their  effect  on  turbulence  fluxes.  The  general  question  we  will  investigate  is  whether  the 
effect  of  these  structures  are  properly  represented  in  an  EDMF  based  PBL  parameterization. 


A  post-doc  (Sandip  Pal)  and  a  graduate  student  (Mark  Sghiatti)  at  UVA  will  contribute  to  the  analysis 
of  the  model  output  and  comparisons  with  observations  during  the  next  year.  Graduate  student  Mark 
Sghiatti  successfully  defended  his  Master  thesis  proposal  entitled  “Spatial  Variability  of  Turbulent 


Kinetic  Energy  and  Turbulent  Fluxes  around  an  Isolated  Mountain”.  He  is  expected  to  graduate  in  fall 
2015. 
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